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ABSTRACT 

Aims. We calculate the energy release associated with a strong first-order phase transition, from normal phase N to an "exotic" 
I , superdense phase S, in a rotating neutron star. Such a phase transition N — >S, accompanied by a density jump p^j — > pg, is 
O ' characterized by Pg/p^ > |(1 + ^b/PnC^), where Pq is the pressure, at which phase transition occurs. Configurations with small 
^ [ S-phase cores are then unstable and collapse into stars with large S-phase cores. The energy release is equal to the difference 
C/3 , in mass-energies between the initial (normal) configuration and the final configuration containing an S-phase core, total stellar 
I ^ I ' baryon mass and angular momentum being kept constant. 

Methods. The calculations of the energy release are based on precise numerical 2-D calculations. Polytropic equations of state 
T-H ■ (EOSs) as well as realistic EOS with strong first-order phase transition due to kaon condensation are used. For polytropic 
^ ' EOSs, a large parameter space is studied. 

Results. For a fixed "overpressure", SP defined as the relative excess of central pressure of collapsing metastable star over the 
pressure of equilibrium first-order phase transition, the energy release Erd does not depend on the stellar angular momentum. 
It coincides with that for nonrotating stars with the same 5P. Therefore, results of 1-D calculations of -Brei {SP) for non-rotating 
stars can be used to predict, with very high precision, the outcome of much harder to perform 2-D calculations for rotating 
f stars with the same SP. This result holds also for 5Pmin < SP < 0, corresponding to phase transitions with climbing over the 
' energy barrier separating metastable N-phase configurations from those with an S-phase core. Such phase transitions could be 
^ realized in the cores of newly born, hot, pulsating neutron stars. 
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1. Introduction a well defined pressure Pq. It is accompanied by a density 

. r . , jump at the phase interface. 

One of the mysteries of neutron stars is the actual struc- 
ture of their superdense cores. Many theories of dense mat- A first-order phase transition in neutron star core is 

ter predict there a phase transition into an "exotic" (i.e., associated with a collapse of initial metastable configura- 

not observed in laboratory) state. Theoretical predictions tions buih of exclusively of non-exotic (normal - N) phase, 

include boson condensation of pions and k aons, and de- into a more compact configuration with a core of the su- 

confinement of quarks (for review sec, e.g., lGlendennin3 perdense (S) exotic phase. At the core edge the pressure is 



[2000: Weber,, . 1999; Haensel et al..„2007,) . Po, and the density undergoes a jump from p„ on the N- 
The first-order phase transitions, accompanied by dis- side to Pg on the S-side. The collapse, called "corequake" , 
continuities in the thermodynamic potentials, seem to be is associated with a release of energy, E.^i. A neutron star 
the most interesting, as far as the the structure and dy- corequake implied by a first-order phase transition in stel- 
namics of neutron stars are concerned. In the simplest l^r core could occur during an evolutionary process, in 
case, one considers states consisting of one pure phase, which central pressure increases. Examples of such pro- 
High degeneracy of the matter constituents implies that cesses are mass accretion and pulsar spin-down. In both 
the effects oftemperature can be neglected. In the thermo- cases, initial and final configurations are rotating. One 
dynamic equilibrium, the phase transition occurs then at assumes that the baryon mass of collapsing star, A'h, is 

conserved (no mass ejection) and that the total angular 

Send offprint requests to: J.L. Zdunik momentum, J, is also conserved (J loss due to radiation 
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during a corequake due to radiation of the electromagnetic 
and gravitational waves is negligible). 

Crucial for the corequake is the value of the param- 
eter A = /?s/p„. If A < Acrit = f(l + -Po//OnC ), then 
the configurations with arbitrarily small cores of the S- 
phase are stable with respect to axisymmetric perturba- 



tions (ISeidov| . ll9Tll;lKa"empfeil . ll98lllHaensel et al.l . ll986t 



Zdunik et al 



19871 ). To get a corequake in an evolving 



neutron star, a metastable core of the N-phase, with 
central pressure Pc > Pq and radius r^^ should form 
first. This core is "overcompressed" , with degree of over- 
compression measured by a dimensionless "overpressure" 
SP = {Pc — Po)/Pq. At some critical value of the over- 
pressure, the S-phase nucleates, and the S-core of radius 
rg forms. For SP — > 0, we have — > 0. This is the case 
of a weak first order phase transition. Up to now, all but 
one numerical calculations were restricted to the spherical 
non-rotating neutron stars (Hacns el fc Proszyhskl. 1982 ; 
Haeiisel et al. . 1986t Zdunik et alT " 1987t Haensel et al 



199d iMuto fc Tatsumil . ll99Cl[ ). An exception is a very re 



cent paper, based on precise 2-D calculations performed 
for w eak first-order phase transitions in rotating neutron 
stars (|Zdunik et all . boOTi ). 



For a strong first-order phase transition (A < Adit) 
in neutron star core, SP — > implies r^ — > r^^niin, 
where rg^min is a sizable fraction of the stellar ra- 
dius. The corequake, accompanying the phase transi- 
tion, is a large-scale phenomenon, with energy release 
Ere\ ~ 10^^ — 10^^ erg. Existing numerical calculations 
were restricte d to spherical non-ro tating neutron stars 
(iMigdal et all. Il979l: iHaensel fc Pr oszvhski. 1 980l. Il982 : 



19821 Il983t iDiaz Alonso 



Kaempfer', '1982"; iBerezin et al 
1983; Bcrczhiani et al., 2003). 

In the present paper we calculate the energy release 
due to a strong first- order phase transition in a rotating 
neutron star. The calculations are done using very precise 
2-D codes and a set of EOSs with strong first-order phase 
transitions. We show, that similarly as in the case of a 



weak first-order phase transition, studied in I Zdunik et al 



()2007( ). the energy released during a corequake depends 
only on the excess of the central pressure of the metastable 
configuration over Pq, and is to a very good approxima- 
tion independent of the angular momentum of collapsing 
star. Moreover, we show that this property holds also for 
corequakes with initial Pc < Pq, which require additional 
energy needed to jump over the energy barrier. 

The paper is organized in the following way. In Sections 
[2] and [3] we introduce notations and we describe general 
properties of the first-order phase transitions in stellar 
core with particular emphasis on the metastability and 
instability of neutron star cores. Analytic models of the 
EOSs with first-order phase transitions, allowing for very 
precise 2-D calculations, are considered in Sect. 14.11 where 
we derive generic properties of the energy release due to 
a first order phase transition at the center of a rotating 
star. In Sect. 14.21 we present our results obtained for a re- 
alistic EOS of normal phase, and we confirm remarkable 
properties of the energy-overpressure relation (i.e., the in- 



dependence from J), obtained in the previous section. In 
Sect.[S]we discuss the decomposition of the energy release 
into kinetic and internal energies. Sect. [6] is devoted to the 
problem of transition of one phase star into stable, two 
phase configuration through the energy barrier. Finally, 
Sect. [7] contains discussion of our results and conclusion. 

2. EOS with a strong first-order phase transition 

At the densities under consideration, all constituents of 
the matter are strongly degenerate, and temperature de- 
pendence of the pressure and energy density can be ne- 
glected. At a given baryon density, nb, energy density of 
the N-phase of matter (including rest energies of parti- 
cles which are matter constituents) is £^{nh) and pressure 
Pjj(nb)- The baryon chemical potential = enthalpy per 
baryon in the N phase, is = (P„ -I- f„)/nb. Similarly, 
one can calculate thermodynamic quantities for the S- 
phase. 





P P P 



Fig. 1. A schematic representation, in the p — P plane, 
of an EOS with a strong first order phase transition, 
Ps/Pn > 1(1 + Po/ Pn^^)- Sohd segments: stable N and 
S phase (in thermodynamic equilibrium). Dash-dot seg- 
ment: metastable N phase. The S phase nucleates at Pnuci, 
which depends on the temperature and the compression 
rate. At Pcrit nucleation of the S phase is instantaneous, 
because the energy barrier, separating the N phase from 
the S phase, vanishes. 



The value of Pq is obtained from the crossing condi- 
tion /^n(P) = /ig(P), which yields also the values of the 
matter densities, p„ and Pg , and the corresponding baryon 
densities, n„ and rig, at the N-S phase coexistence inter- 
face. These parameters are obtained assuming thermody- 
namic equilibrium. A schematic plot of EOS of matter 
with a strong first order phase transition N — >S, in the 
vicinity of the phase transition point, is plotted in Fig. 
[H Solid segment of the N-phase corresponds to the sta- 
ble N-phase states. For pressures above Pq, the N phase 
becomes metastable with respect to the conversion into 
the S phase. The S-phase can appear through the nucle- 
ation process - a spontaneous formation of the S-phase 
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droplets. However, an energy barrier resulting from the 
surface tension at the N-S interface delays the nucleation 
for a time identified with a lifetime of the metastable state 
''nuci- The value of t^ucI decreases sharply with P > Pq, 
and drops to zero at some Pcritj where the energy bar- 
rier separating the S-state from the N-one vanishes. For 
P > Pcrit the N phase is simply unstable and converts 
with no delay into the S phase. Calculations of Tnuci in 
dense neutron star core metastable with respect to the 
pion c onde nsation was performed by iHaensel fc Schaeffer 



pion c onae nsation was pertor r ned Dy itiaensei &: acnaeitei 
i|l982') and i Muto fc Tatsumil (|l99(X ). The case of nucle 
ation of quark matter was studied by llida fc Satol ()l997l . 

while nucleation of the kaon condensate was con- 
sidered byfN orseni (j2002l ). Discussion of Tnuci in an accret- 



ing or spinni ng down neutron sta r with increasing P^ was 
presented by IZdunik et all (|2007l ); a spe cific case of a hot 
accret ing neutron star was considered bv lBerezhiani et al 
( 2003h . 



3. Families of neutron stars 

For a given EOS, nonrotating hydrostatic equilibrium con- 
figurations of neutron stars form a one-parameter family, 
the parameter being, e.g., central pressure Pc- Notice, that 
Pc is preferred over pc, because pressure is strictly contin- 
uous and monotonous in stellar interior, while density can 
suffer discontinuities. For an EOS with a first order phase 
transition, and without allowing for metastability of the 
N phase, neutron stars form two families: a family of stars 
composed solely of the N-phase {C}, and a family of those 
having an S-phase core {C*}. 

It is known since longtime that configuration {C*} with 
a small S- phase c o re are unstable for A > Aerit = 1(1 + 
(ISeidov 



19861: 



Zdunik et al 



19711: iKaempferi . 119811 : IHaensel et al 



1983). Therefore, the topology of the 



complete set of stable hydrostatic configurations of non- 
rotating neutron stars, parametrized by Pc, and plotted, 
e.g., in the mass-radius plane, depends on the value of 
A. Namely, for A < Acrit, the sum of {C} and {C*} is 
continuous, while for A > Acrit it is not (i.e., the families 
{C} and {C*} are disjoint). Recently, it has been shown, 
that this property is generic for an EOS with a phase 
transition, and holds also for rigidly rotating neutron stars 



( Zdunik et all . l2006l ) 



An important global parameter for a hydrostatic equi- 
librium configurations is its baryon mass, il/b. It is de- 
fined as the baryon number of the star. Ah, multiplied by 
a "mass of a baryon" , toq, defined as the 1/56 of the mass 
of the ^^Fe atom: toq — 1.6586 x 10~^^ g. During evolution 
of an isolated neutron star, including the phase transitions 
in its interior, Mb remains strictly constant. 

We define the equatorial radius of an axisymmetric 
neutron star as the proper length of the equator divided by 
2tt. Examples of M^-Rcq curves for non-rotating and ro- 
tating neutron stars without a phase transition {C}, with 
a weak first-order phase transition {0^^,^^.}, and a strong 
first-order phase transition {CJ^j-cig}, are shown in Fig. 
[21 Dotted segments correspond to unstable configurations 




Fig. 2. Baryon mass vs. equatorial radius for hydrostatic 
equilibrium configurations calculated for three types of 
EOS of dense matter, described in the text. Solid line - 
stable; dotted line - unstable configurations. Thick lines 
correspond to the non-rotating models, thin lines to the 
rigidly rotating configurations with a fixed total stellar 
angular momentum J = 1.2 x GMq/c. 



(instability with respect to the axisymmetric perturba- 
tions). As one sees, fast rotation significantly changes the 
Afb(Pcq) dependence, e.g., by increasing the mass and the 
radius of the configuration with Pc = Pq denoted by Co . 

General relations between models calculated for three 
types of the EOS can be formulated. At a given Mb, 
■Req(Q„ng) < ^oq(C*cak) < -Rcq(C). Morcover, maxi- 
mum allowable baryon masses satisfy Mb,max(C*tronK) < 



M 



b,max(C*ga]^) < Mb,max(C); the Same inequalities are 
valid for maximum allowable gravitational mass, Mmax- 

Differences in the mass-radius behavior are most pro- 
nounced in the vicinity of configuration Cq, with Pc = Pq- 
This region, bounded by a rectangle in Fig. [2l is shown 
in Fig. [21 where the arrows connect configuration with 
same Mb. For simplicity, we first consider non-rotating 
stars. For A < Adit, configurations with Pc > Pq form a 
monotonous branch C*cak- -f'o^ > Acrit, a segment with 
Pq < Pc < Pc*i„in consists of Configurations C*^_^^^^ with 
dM*/d/9c < 0, which are therefore unstable with respect 
to radial perturbations. Therefore, for A > Acrit the stable 
branches {C} and {C*} are disjoint. In both cases ("weak" 
and "strong"), central density jumps from p„ to when 
passing from the C branch to the C* branch. All these 
properties were derived long time ago for non-rotating 
neutron stars. Recently it has been shown that they hold 
also for phase transitions in rigidly rotating neutrons stars 
when the stationary axially symmetric families C and C* 
contain configuration of a fixed stellar angular momen- 
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turn J ( Zdunik et al. . 2006f) . The standard non- rotating 
case corresponds to J = 0. 



I J, f 




Fig. 3. Zoomed fragment of Fig. [2l in the vicinity of the 
phase transition. For other explanations see the text. 



4. Calculation of the energy release 

We restrict ourselves to axially symmetric, rigidly rotating 
neutron stars in hydrostatic equilibrium. In what follows 
by "radius" we mean the equatorial circumferential radius, 

We assume that at a central pressure Pc = -Pnuci the 
nucleation of the S phase in an overcompressed core, of 
a small radius r-^ (of configuration C) initiates a strong 
first-order phase transition. This leads to formation of a 
large S-phase core of radius in a new configuration C*, 
as shown in Fig. [4] 

We compare the hydrostatic equilibria of neutron 
stars corresponding to the EOSs with and without the 
phase transition using the numerical library LORENE 
(Ihttp : //www . lorene . obspm . f r) , obtaining the axisym- 
metric , ri gidly rotati ng so lutions of Einstein equations 



as in IZdunik et all (|2007t) 



The accuracy of the solu- 
tion, measured with the general re lativistic virial theorem 
( Bonazzola fc Gourgoulhon . 1994[ l is typically 10~^. 

The neutron-star models can be labeled by the cen- 
tral pressure Pc (central density is not continuous!) and 
rotational frequency / — r2/27r. These parameters are nat- 
ural from the point of view of numerical calculations. But 
we can introduce another parametrization, more useful for 
other purposes. In order to study the stability of rotating 
stars, a better choice is the central pressure, Pc, and the 
total angular momentum of the star, J. 

We additionally assume that the transition of the star 
from a one-phase configuration to the configuration with a 
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Fig. 4. Transition from a one-phase configuration C with 
a meta-stablc core of radius to a two-phase configu- 
ration C* with S-phase core with a radius . These two 
configurations have the same baryon mass Afb = and 
total angular momentum J = J* . 



dense core of the S-phase takes place at fixed baryon mass 
A/b (no matter ejection) and fixed total angular momen- 
tum of the star J (loss of J due to the electromagnetic 
or gravitational radiation is neglected). The energy re- 
lease during transition C(Mb,/) — > C*{Mb,f*) is there- 
fore calculated from the change of the stellar mass-energy 
during this process. 



E,,, = c2 [M{C) - M{C*)] 



Mi, J 



(1) 



The precision of the determination of Ej-d depends on 
the numerical accuracy of calculations of the configura- 
tions with the same J and Mb. Taking into account that 
Proi is 3 4 orders of magnitude smaller than M (for the 
overpressures ranging from a few to 20 percent) this pre- 
cision is typically better than 1%. 

4.1. Energy release for polytropic EOSs 

The use of the polytropic EOSs for the N and S phases 
not only guarantees very high precision of numerical cal- 
culation, but opens also a possibility of the exploration 
of wide region of the parameter space. Description of 
the polytropic EOSs and their application to relativis- 
tic stars with phase transitions was prese nted in detail 



in our previous publicati ons in this series (jBeiger et al 



2005 : Zdunik et al. . 20061 ). In Fig.[6]we present the energy 



release as a function of , for several values of the angular 
momentum of the metastable configurations C(J). As we 
see the energy release corresponding to a given value of 
Tg depends rather strongly on the total angular momen- 
tum. For example, the value of E^d at = 6.4 km, for 
J — (nonrotating star), is twice that for J = 1.0 GMq/c. 
Moreover, let us notice that the value of E^c\ for minimum 
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Fig. 5. Total baryon mass Mb of hydrostatic stellar config- 
uration, versus central pressure Pc, at fixed stellar angular 
momentum J, for the EOS depicted in Fig. [T] Solid line 
denotes stable states, dash-dot line - the states which are 
metastable with respect to the N — >S transition. For a 
central pressure Pnuci the S-phase nucleates in the super- 
compressed core of configuration C, and this results in a 
transition C — > C* into a stable configuration with a S- 
phase core and central pressure P* . Both configurations C 
and C* have the same baryon mass M\y. 

Table 1. The parameters of the polytropic EOSs. The 
basic EOS for the N phase, named PolN, is P = 
Kj^c'^{nh/ni)'^r^ , where ni = O.f fm^^ and = 
4.15 fO"'^^ erg cm~'^, and 7„ = 2.5. Construction of the 
EOS for the S-phase and of the first order ph ase transition 
at P = Pq is based on the Appendix A of IZdunik et al. 

The transition point is the same for all three 
EOSs, and is defined by Pq = 3.686 lO^Vg cm'^ = 
0.25 fm~^, and = 4.4210^'*g cm"^. The density jump 
corresponding to the phase transition is defined by the 
parameters A — Ps/Pti and An = n^/n^, connected by 
the relation A = 1 + (An - 1)(1 + Po/p^c^)- The S-phase 
EOS is a polytrope with the parameters and 7s . The 
non-rotating reference configuration Cq has M = Mq = 
1.38 Mq and R = Ro ~ 15.2 km. Phase transition in 
PolWl is weak (A < Adit), while those in PolSl and PolS2 
are strong (A > Acrit)- 



EOS 


(erg cm 


7s 

(erg cm ^) 


A 


An 


PolW 


5.11 10'' 


3.5 


1.44 


1.4 


PolSl 


1.602 W 


4 


1.66 


1.6 


PolS2 


1.257 W 


4 


1.77 


1.7 



Tg (at fixed J) is (nearly) independent of J. This refiects 
independence of E-^f,i{5P = 0) from J. 

In Fig. [7] we present the energy release as a function 
of the overpressure of the metastable N phase in the cen- 
ter of the metastable star C(Mb, J), for several values of 
J. As we already stressed, the value of Pnuci (or 5P) can 
be determined from microscopic considerations, combined 
with physical conditions prevailing at the star center as 
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Fig. 6. (Color online) The energy release due to the core- 
quake of rotating neutron star as a function of the equa- 
torial radius of the S-phase core, . Calculations are 
performed for the EOS PolSl from Table 1. Different 
curves correspond to the different values of total angu- 
lar momentum of rotating star, fixed along each curve, 
J = J GM%lc = (0,0.1, . . . , 1.3) X GMl/c, from the 
left-most curve to the right-most curve. 



well as with the time evolution rate. Having determined 
Pnuci, we can determine the energy release, E'rci, due to the 
corequake C(Mb, J) — > C*(Afb, J), where the metastable 
one-phase configuration, and the final two-phase configu- 
ration, have the same values of the baryon mass Afb and 
total angular momentum J. 

As we see in Fig. [71 the energy release in a collapse of a 
rotating star is independent from the angular momentum 
of collapsing configuration, and depends exclusively on the 
degree of metastability of the N phase at the stellar cen- 
ter (departure of matter from chemical equilibrium) , mea- 
sured by the overpressure 5P. Consequently, to obtain the 
energy release associated with a corequake of a rotating 
neutron star, it is sufficient to know the value of Erd for 
a non-rotating star of the same central overpressure. It 
should be stressed that the configurations C(Mb, J) and 
C*(Mb, J), considered in this section, are the fast rotat- 
ing ones, those with largest J are close to the Keplerian 
(mass shedding) limit. This is visualized in Fig. [H where 
we plotted the oblateness of the star and the kinetic to 
potential energy ratio. And still, in spite of fast rotation 
and large oblateness, the energy release is the same as in a 
non-rotating star of the same initial central overpressure. 

Summarizin g, a remarkab l e inde pendence of E-^d from 
J, obtained in IZdunik et al. I (|2006h for A < Acrit, when 
the small-core approximations were valid, holds also for 
A > Acrit, where perturbative arguments cannot be used. 
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Fig. 7. (Color online) The energy release due to the mini- 
collapse of rotating neutron star as a function of the over- 
pressure 5P of the N phase of the matter in the center 
of the star, for three EOSs from Table 1. The bottom 
curve is calculated for the PolW EOS wit a weak first- 
order phase transition, so that -Erci(O) — 0. The middle 
curve and the top curve are for the PolSl and PolS2 
EOSs, respectively, with strong first-order phase transi- 
tions. The points of different color correspond to the dif- 
ferent values of total angular momentum of rotating star, 
J ^ J GMl/c = (0, 0.1, . . . , 1.3) X GMl/c. For a given 
EOS, results for all rotating configurations can be very 
well approximated by a single curve, independent of J. 



4.2. Realistic EOS - an example 

In the present section we check the general validity 
of results obtained for polytropic EOS and discussed 
in Sect. 14.11 by performing numerical calculations for 
realistic EOSs with a strong first-order phase transi- 
tion. For the EOS of the crust we took the model of 

The constituents of the N 



Douchin fc Haensell (|200lh . 
phase of the core were neutrons, protons, electrons, and 
muons. Nucleon component was described using the rela- 
tivistic mea n-field model with scala r self-c oupling, con- 
structed by Zimanvi fc Moszkowski ( 1990l ). The values 
of the meson- nucleon coupling constants were go-/ ma- = 
3.122 fm, = 2.1954 fm, gp/rup = 2.1888 fm. 

The dimensionless coefficients in the cubic and quartic 
terms in scalar self-coupling were b = —6.418 x lO^'^ 
and c = 2.968 x 10"'^, respectively. As an example of 
the S-phase we considered the kaon-condensed matter. We 
constructed a specific dense matter model with a strong 
first-order N — >S transition implied by kaon condensa- 
tion. Cou pling of kaons to nucleons is des c ribed by the 
model of iGlendenning &: Schaffner- BiehchI (Il999l) .' with 
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Fig. 8. Top panel: the ratio of polar radial coordinate to 
the equatorial radial coordinate ratio. Bottom panel: the 
ratio of the kinetic energy, T and the absolute value of 
the potential energy, W, for the reference stellar configu- 
rations (central pressure Pq, Table 1) consisting of the N 
phase of dense matter, described by the polytropic EOS, 
PolN, of Table 1. Large dots correspond to the values of 
the total stellar angular momentum, J ~ JGMq/c = 
(0.1, . . . , 1.5) X GM^/c, which were used in Figs. M El 



[/li." = -115 MeV. Resulting EOS is shown in Fig. H 
The phase transition is the strong one, with A > Acrit (see 
the caption to Fig. [5]) . 

In Fig. [TO] we show the energy release due to the 
C{Mh,J) — > C*(Mb, J) transition, versus overpressure. 
The values obtained for different values of J are marked 
with different color and symbols. To an even better ap- 
proximation than for the polytropic models, all color 
points lie along the same line. For a given overpressure 
5P, the energy release does not depend on J of collapsing 
metastable configuration. This property holds for a broad 
range of of stellar angular momentum, J — (0.1, . . . , 1.0) x 
GM|/c. 



5. Decomposition of the energy release 

As was already said, we restrict ourselves to the sta- 
tionary, axisymmetric states of rotating neutron stars. In 
the Newtonian theory energy of a rigidly rotating axially 
symmetric star (body) is easily decomposed into kinetic 
energy of rotation, T, and internal energy, U, which is 
the the total energy of the star measured in the stel- 
lar (body) reference system. Kinetic energy of rotation is 
T — 5 Jfi ~ 5^^^' where the moment of inertia / = J/fl. 
On the other hand, the total energy measured in "the lab- 
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Fig. 9. EOS with first-order phase transition, due to 
kaon , described and used in Sect. 14.21 in the present 
paper. The phase transition occurs at Pq = 4.604 x 
lO^'^ erg cm-3. On the N-phase side, n„ = 0.3375 fm-^ 
and — 5.973 x lO^'' g cm^'^. On the kaon-condensed S- 
side, = 0.6122 fm^^ g^^j ^ i^j25 x 10^^ g cm'^. 
Therefore, A — p^/p^ = 1-884, which is greater than 
Acrit = |(1 + Po/Pnc') - 1.629. 



oratory frame" , E, is related to U hy U = E — Jil. For a 
rigid rotation U = E — ^J^- 

In general relativity, all kinds of energies sum up to 
give the stellar gravitational mass, M — E/c^, which is the 
source of the space-time curvature. Therefore, the decom- 
position of Mc^ into T and U is ambiguous. Here, we will 
use a standard Newtonian-hke formula, T = ■^J fl, where 



both J and are well defined quantities (jFriedman et al 
198i). 



The transition C{Mh,J) — > C*{Mh,J) is accompa- 
nied by a spin-up of neutron star, and an increase of its 
kinetic energy. Our results show, that while i^rei is to a 
very good approximation independent from J, the kinetic 
energy increase AT = {T* - T)m^,j = ^Ji^* - n)M^,j 
grows rapidly with J. This is visuahzed in Fig. [TOl The 
contribution to the energy release resulting from the de- 
crease of the internal energy of the star will be denoted 
by eI'^*'^ . It is given by 



E. 



(int) 
rol 



{SP, J) - -Jin* - + E,,,{SP) 



Therefore, the energy release in the star reference system 
is greater than the total energy release, measured by a 
distant observer. The difference increases rapidly with J, 
see Fig. [TUl 



^ 10 - 



< 




Fig. 10. (Color online) The total energy release, E^rci, due 
to the collapse of a rotating neutron star, implied by kaon 
condensation, as a function of the metastabihty (overpres- 
sure) of the normal phase of the matter in the center of the 
star for the EOS in Fig.[9l The points of different color cor- 
respond to the different values of total angular momentum 
of rotating star. The results for all rotating configurations 
can be very well approximated by one curve. Solid lines: 
kinetic energy increase, AT — T* ~ T, for fixed values 
of J = JGMl/c = (0.1, . . . , 0.8) x GMl/c. Bottom hue 
is for J = 0.1, and top line for J = 0.8. Note the dif- 



ferent definition of the sign of T'rci = (M 
AT = T* - T. 



M*)c^ and 



6. Phase transition with climbing over the energy 
barrier 



In the case of a strong first-order phase transition, config- 
urations with Tcrit > > ^0 are not the only ones which 
are metastable with respect to the N — >S transition. In 
order to look for additional metastable N-phase configu- 
rations we plotted, in Fig. 1111 the vicinity of the reference 
configuration Co in the Afb — R plane: for Mb (Co) > Mb > 
Mb(C^i„) we have three equilibrium configurations with 
a given Mb. Consider a triplet of equilibrium configura- 
tions CiCi'C*. Two of them, Ci and Cjf, are stable (local 
minimum of energy), and one C^' is unstable (local max- 
imum of energy). Notice that the N-phase configurations 
on the CoCmin segment are characterized by negative over- 
(2) pressure (which might then be called "underpressure") 

SP = Pnucl/Po - K 0. 

In Fig.[TT]we indicated, by horizontal lines, several ex- 
amples of transitions between equilibrium configurations 
with the same Mb. To understand the nature of these 
transitions, we plotted in Fig. [12] the energy release, Eq. 
([1]), associated with a transition between a pair of con- 
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Fig. 11. Three families of neutron stars in the Mb — R, 
plane and the transition Ci 
energy barrier. 



oq 



Ci via jumping over the 



figuration, versus the negative overpressure. We see that 
the functional dependence Erci{SP), together with (quite 
precise) independence from J, continue smoothly into the 
region of negative SP. Notice that point A in Fig. [12] cor- 
responds to Co — > Cq. As far as the transitions Ci — > 
(segment BD in Fig. [T^ are concerned, they are always 
associated with E^^i < 0, i.e. to make them the star should 
gain (absorb) energy instead of releasing it. The necessary 
energy input in B — ji^rcil, and it reaches a maximum at 
point D, corresponding to Cmin — * ^min- Therefore, in 
order to get to C*' by forming a small S-phase core, the 
system has to climb over the energy barrier. Then, config- 
uration Ci (which is unstable) can collapse into with a 
large core, and this collapse is associated with an energy 
release. In this way the star reaches the global minimum 
of M at fixed Mb. However, this is the case provided the 
transition takes place above the horizontal line C2C2. 

Summarizing, N-phase configurations on the C0C2 seg- 
ment are metastable with respect to transition to large S- 
phase core configuration on the CgCj segment. However, 
such a transition requires climbing over the energy bar- 
rier of a height B. Let consider a specific example with 
5P = -0.02. Using Fig. [El we see that B ~ 10"*^ erg while 
i?rei = 3 X 10^*^ erg. Generally, for this model of phase tran- 
sition we have i^rci S> B for underpressures 5P > —0.03. 
The excitation energy, i?oxc, contained in radial pulsa- 
tions, scales as the square of relative amplitude 6R/R. 
For the fundamental mode, Ecxc ~ 10^^ {5R/ R)"^ erg. 
Therefore, E^xc exceeds 10^^ erg as soon as SR/R > 0.01, 
a condition that is easy to satisfy by an newly born neu- 
tron star. 

However, a second obstacle for the collapse of configu- 
ration with a negative SP should be pointed out. Namely, 



Fig. 12. (Color online) Energy release, Eq. ([T]), versus 
overpressure SP — {Pc — Po)/Pq, for transitions Ci — > 
Ci and Ci — > CJ'' visualized in the Mb — R plane in 
Fig. [11] Points of different color correspond to transi- 
tions with different angular momentum J = JGMq/c^ = 
(0,0.1, . . . ,0.9) X GMl/c'^. Solid line - results for J = 
(non-rotating stars). Point A corresponds to Co — > Cq. 
Point B corresponds to Co — > Cq' = Co, point D to 
Cmin — * C,* i„, and point C to C2 — > C^. 



apart from the energy condition allowing climbing over an 
energy barrier, there is a timescale condition: there should 
be enough time to form the S-phase core. This latter con- 
dition may be more difficult to fulfill than the former one, 
particularly if there is a need to create strangeness, like 
in the kaon condensation or in the formation of three- 
fiavor u-d-s quark matter from a deconfined two-fiavor (u- 
d) state. Once again, favorable conditions for Ci — > C{ 
with climbing over Ci' could exist in a newborn neutron 
star. A neutron star born in gravitational collapse not only 
pulsates, with pulsational energy much greater than B, but 
additionally a high temperature ~ 10^^ K in the stellar 
core can allow for a rapid nucleation of the S-phase. 

7. Discussion and conclusions 

The most important result of the present paper is that 
the total energy release, associated with a strong first- 
order phase transition at the center of a rotating neutron 
star, does depend only on the overpressure at the center of 
the metastable configuration and is independent from the 
star rotation rate. This result holds even for fast stellar 
rotation, when the star shape deviates significantly from 
sphericity, and for overpressures as high as (10-20)%. This 
property is of great practical importance. Namely, it im- 
plies that the calculation of the energy release for a given 
overpressure, requiring very high precision 2-D calcula- 
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tions to guarantee Mb = M^, can be replaced by a sim- 
ple calculation of non-rotating spherical stars. The inde- 
pendence of the energy release on the rotation should be 
treated as a result of numeric;al c;alculations and is subject 
to the numerical accuracy of the stellar parameters deter- 
mination. Strictly speaking our numerical results indicate 
that, if the energy release depends on the rotation, this 
dependence is extremely weak. Namely, maximum devia- 
tion from the nonrotating value for a given overpressure is 
of the order of 1%, which actually is the precision of our 
numerical calculations. 

We studied stability of configurations with central 
pressure below that for the equilibrium phase transition. 
If the initial state of neutron star is excited, e.g., is pul- 
sating, then the formation of a large dense phase core is 
possible, but it requires climbing over the energy barrier 
associated with formation of a small core. Excitation en- 
ergy has to be larger than the height of the energy bar- 
rier. Additionally, if a phase transition is connected with 
a change of strangeness per baryon, then the temperature 
has to be high enough to make strangeness production 
sufficiently rapid. Such conditions might be realized in the 
cores of newly born neutron stars. 

Note that the energy release Ere\ ~ 10^^ — 10®^ erg, is 
an absolute upper bound on the energy which can released 
as a result of a phase transition at the star center. The 
energy AE can be shared between, e.g., stellar pulsations, 
gravitational radiation, heating of stellar interior. X-ray 
emission from neutron star surface, and even a gamma- 
ray burst. 
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